The enumeration of simple permutations 
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Abstract 

A simple permutation is one which maps no proper non-singleton 
interval onto an interval. We consider the enumeration of simple per- 
mutations from several aspects. Our results include a straightforward 
relationship between the ordinary generating function for simple per- 
mutations and that for all permutations, that the coefficients of this 
series are not P-recursive, an asymptotic expansion for these coeffi- 
cients, and a number of congruence results. 
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1 Introduction and definitions 

The permutation 2647513 maps the interval 2. .5 onto the interval 4. .7. In 
other words, it has a segment (set of consecutive positions) whose values 
form a range (set of consecutive values). Such a segment is called a block of 
the permutation. Every permutation has singleton blocks, together with the 
block l..n. If these are the only blocks the permutation is called simple. For 
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example, 58317462 is simple and the simple permutations of length up to 5 
are as follows:. 



Length Simple permutations 

1 1 

2 12, 21 

3 None 

4 2413, 3142 

5 24153, 25314, 31524, 35142, 41352, 42513 

Simple permutations have recently had important applications in the study 
of pattern closed classes of permutations p. 

Let Sn denote the number of simple permutations of length n. We shall 
be concerned with properties of the sequence Consider the ordinary 

generating functions: 

oo 

k=l 

oo 

S{x) = Y.^kX^- 

fc=4 

We start S{x) from x^ because simple permutations of length 1 and 2 need 
special treatment. Later in this section we will see that the coefficients of 5* 
differ from those of — (functional inverse, not reciprocal) alternately by 
2 and —2. The coefficients of were considered by Comtet j3 p. 171] 

without any combinatorial interpretation. The sequence of absolute values 
of these coefficients appears as sequence A059372 of [12], and the first few 
terms are: 

1, 2, 2, 4, 4, 48, 336, 2928, 28144, 298528, 3454432, 43286528. 

So we shall see that the numbers s„ are: 

1, 2, 0, 2, 6, 46, 338, 2926, 28146, 298526, 3454434, 43286526. 

In section 12] we shall prove that (s„) is not P-recursive (it cannot be defined 
by a linear recurrence with polynomial coefficients). In section JS] we derive 
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Figure 1: A block decomposition of 67183524. The pattern of the block 
decomposition is the permutation whose graph is defined by the occupied 
cells, namely 3142. Within each occupied cell, the individual blocks also 
define permutations namely 12, 1, 1, and 2413. 

the asymptotic behaviour of s„ (the main term is nl/e^) and section |3] gives 
various congruences satisfied by the numbers s„. 

In the remainder of this section we derive a structure theorem that shows 
how arbitrary permutations are built from simple ones, and read off from it 
equations satisfied by generating functions. We begin with some terminology 
and notation that will be used throughout. 

A block decomposition of a permutation a is a partition of a into blocks. Of 
course, if a is simple there will only be the two trivial block decompositions. 
An example of a non-trivial decomposition is a = 67183524 with blocks 
(67)(1)(8)(3524). 

Given a block decomposition of a, its pattern is the permutation defined by 
the relative order of the blocks. In the example above, the pattern of the block 
decomposition (67) (1) (8) (3524) is 3142. We may think of the permutation 
67183524 as being constructed from the permutation 3142 by inflating each 
of the elements into a block, in this case the blocks 12, 1, 1, and 2413 (we 
view each block as a permutation in its own right). We write: 

67183524 = (3142) [12, 1, 1, 2413]. 

This example is further illustrated in Figure ^ The inflation procedure is an 
instance of the wreath product for permutations 

A permutation which cannot be written in the form (12) [a, /5] is called plus 
indecomposable, and one which cannot be written in the form (21) [a,/?] is 
called minus indecomposable. Let in denote the number of plus indecompos- 
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able permutations of length n. The number of minus indecomposable per- 
mutations of length n is also i„ as is easily seen by considering the bijcction 
on permutations of length n which sends tt to tt' where 7r'(t) = n + 1 — 7r(t). 

Theorem 1 For every non-singleton permutation n there exists a unique 
simple non-singleton permutation a, and permutations cki, 0:2, . . . , ccfc such 
that 

Moreover, if cr ^ 12, 21 then cti, a2, . . . ,ak are also uniquely determined. If 
a = 12 (respectively 21) then ai and (i-TC uniquely determined subject to 
the additional condition that ai be plus (respectively minus) indecomposable. 

The caveat added for the case where a = 12 (or 21) is necessary, as is easily 
seen by considering tt = 123. This can be decomposed as (12) [1, 12] or as 
(12) [12, 1]. However, only the former decomposition has a plus indecompos- 
able first part. 

Proof: We first of all suppose that tt has two distinct maximal proper 
blocks A and B that have a non-empty intersection. Then, as the union of 
intersecting segments is a segment and the union of intersecting ranges is a 
range, ^4 U S is a block. Because of the maximality, AU B = [n]. But it is 
also clear that A cannot be an interior segment of [n] nor can it define an 
interior range. In other words we have 

TT = a[a, /?] 

where o" = 12 or o" = 21. These two possibilities are obviously mutually 
exclusive. In either case consider all decompositions of tt as o'[7,5]. The 
intersection of their 7 parts is also the 7 part of a decomposition of this type. 
So there is a unique such decomposition with smallest 7 part. Clearly, this 
part is plus indecomposable in the case o" = 12 and minus indecomposable if 
(7 = 21. 

We next suppose that every pair of distinct maximal blocks has empty inter- 
section. Obviously, then the maximal blocks form a block decomposition of 
TT and this decomposition must be coarser than every other block decompo- 
sition of TT. It follows that this decomposition is the only one whose pattern 
(7 is simple and so we obtain the unique representation claimed for tt. ■ 
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We shall shortly see that this theorem gives relations between the following 
three generating functions: 



oo 




k\x' 



.k. 



k=l 



oo 




.k. 



oo 



S{x) 



fc=4 



Note that our generating functions are all taken to have zero constant term. 
This slightly unconventional choice turns out to be algebraically convenient 
at several points. 

From Theorem ^ it is easy to see that there is a one to one correspondence 
between the collection of all permutations with length at least 2 and the 
collection of sequences: 



Here a may be any simple permutation of length k > 2, and if a 7^ 12, 21 
then ai through are arbitrary permutations, while if a = 12 (respectively 
21), ai is plus- indecomposable (respectively minus indecomposable) and ^2 
is arbitrary. 

This correspondence, together with the earlier observation that the numbers 
of plus and minus indecomposable permutations of length n are the same, 
translates naturally into the following equation: 



However, since a plus indecomposable permutation cannot correspond to 
a sequence beginning with 12, while all other sequences do represent plus 
indecomposables, it is also clear from the correspondence that 



Solving this latter equation for /, and then substituting in equation (H)) before 
solving for S o F gives: 



(a, ai, ^2, . . . , ak). 




(1) 



/(x) = x + I{x)F{x) + {So F){x). 
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Now letting t = F{x) we obtain: 

Sit)=t-^^-F^-'\t). (2) 

We can also obtain an equation for the ordinary generating function of plus 
indecomposable permutations through the observation that every permuta- 
tion decomposes into a sequence of plus indecomposable permutations so 

. I(x) 
Fix) = 



1 - I{x) 
or equivalently 

Denoting the coefficient of in F^^^\t) by Com„ (in reference to Comtet 
who initiated the consideration of this sequence) we obtain directly from 
equation (j^J the simple relationship that for n > 4: 

s„ = -Com„ + (-1)"+^ ■ 2. 



2 Non P-recursiveness 



A sequence of numbers (a„) is called P-recursive if it satisfies a linear re- 
currence with polynomial coefficients. A power series is called D-finite if 
it satisfies a linear differential equation with polynomial coefficients. A 
sequence (a„) is P-recursive if and only if its ordinary generating func- 
tion A{x) = J2n 

D-finite. More information on D-finiteness and 
P-recursiveness can be found in Stanley [IHl Chapter 6]. If = n! then 
ttn — nun-i = 0, and thus the sequence (n!) is P-recursive. We show that on 
the other hand neither sequence (in) nor (s„) is P-recursive. By (j21), instead 
of the latter sequence we can work with (Com„). 

Proposition 2 The power series I (x) andC{x) = F^^^'^(x) = ComfcX*^ 
satisfy the differential equations 

I' = ~x'^I^ + {x-^ + x-^)I-x-^; 

c = 

x-{l +x)C' 



6 



Proof: It follows from the recurrence for n! that F{x) satisfies x + xF + 
x'^F' = F. Thus F' = Combining this with F = J/ (1 - /) 

we obtain the differential equation for I{x). Similarly, C = 1/F'(C) = 
C"^ I ((1 — C)x — C) which is the differential equation for C(x). ■ 

Klazar |HI used the following method to show that a sequence (a„) is not P- 
recursive. Suppose that the ordinary generating function A[x) is non-analytic 
and satisfies a first order differential equation A! = R{x, A) where R is some 
expression. Differentiating this relationship and replacing A' by R{x,A), 
the derivatives of A are expressed as A^''^ = Rk{x,A); Rq{x,A) = A and 
Ri{x, A) = R{x, A). Substituting Rk{x, A) in the equation of D-finiteness 

boA + biA' + h^A" + ■■■ + hsA^'^ = 0, 

where s > 1, 6j G C(a;) and hs 7^ 0, we get a non-differential equation 
Hk=Q bkRk{x, A) = 0. If -R is such that the expressions Rq, Ri, R2, ■ ■ ■ are (i) 
analytic or even algebraic and (ii) linearly independent over C(x), we have 
a nontrivial analytic equation for A. This implies that A is analytic (see 
Klazar's paper 'E] for more details) which is a contradiction. So A cannot be 
D-finite and the sequence of its coefficients cannot be P-recursive. 

To state the result of jH] precisely, we remind the reader that a power series 
R{x,y) e C[[x, I/]] is analytic if it absolutely converges in a neighborhood of 
the origin and that R{x,y) G C{{x,y)) is an analytic Laurent series if, for 
some positive integer k, {xy)^R{x,y) G C[[x, y]] is analytic. Theorem 1 of 
[S] says that if A G C[[x]] is non-analytic, R{x,y) G C((x, y)) is analytic, 
A' = R{x,A), and R contains at least one monomial ax^y^ , a 7^ 0, with 
j < 0, then A is not D-finite. This result applies directly neither to I{x) nor 
C{x) (see Proposition 121) because in the case of /(x) the last condition on R 
is not satisfied and in the case of C{x) the right hand side R even cannot be 
expanded ClS db Laurent series. 

However, the substitution x — (1 + x)C{x) = 6{x) transforms the second 
differential equation of Proposition |21 into 

x^ 1 l + 2x 

" — ■ 7T + ~; • 

1 + X 9 1 + X 

Now all conditions are satisfied {F{x) is clearly non-analytic which implies 
that C(x) and ^(x) are non- analytic) and thus ^(x) is not D-finite by Theorem 
1 of 1^. The dependence of C(x) and S{x) on 9{x) and the fact that D-finite 
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power series form a C(x)-algebra (^21 Theorem 6.4.9]) shows that neither 
C{x) nor S{x) is D-finite. 

In order to deal with the case of I{x), we use this opportunity to complement 
Theorem 1 of jH] in which R & C{{x,y)) hj the following theorem which treats 
the case R G C{x,y). Neither of the theorems subsumes the other because 
not every rational function in x and y can be represented by an element of 
C{{x,y)) (as we have seen) and, of course, not every Laurent series sums up 
to a rational function. However, the next theorem seems to be more useful 
because in both examples in |H] and both examples here the right hand side 
R{x,y) is, in fact, a rational function. 

Theorem 3 Let P,Q E C[x,y] he two nonzero coprime polynomials and 
A G C[[a;]] be a non-analytic power series which satisfies the differential 
equation 

Q{x,Ay 

IfdegyQ = and degyP < 1 then A is, trivially, D-finite. In all remaining 
cases A is not D-finite. 

Proof: The first claim is clear. If deg^^ Q = and r = deg^^ P > 2 then 
A' = aQ-\-aiA-\-- ■ ■ + arA^ where aj G C(a;), r > 2, and 7^ 0. Differentiation 
by X gives 

A^^^ = Rk{x, A) = ao,k + ai^kA + ■■■ + akr-k+i,kA''^~'"^^ 

where Oj j G C(x) and 

akr-k+i,k = r{2r - l)(3r - 2) . . . ((A; - l)r - A; + 2)a,' ^ 0. 

Thus Rk{x, y) G C(a;) [y] have y-degrees kr — k + 1, A; = 0, 1, 2, . . ., which is 
for r > 2 a strictly increasing sequence. Therefore Rq, Ri, R2, ■ ■ ■ are linearly 
independent over C(x) and, by the above discussion, A is not D-finite. 

In the remaining case deg^Q > 1. Differentiation of A' = R{x,A) = 
P{x, A)/Q{x, A) by x gives A^'^'^ = Rk{x,A) where Rk{x,y) G C{x,y). For 
example, 

(P,. + PyRl)Q - P{Q, + QyRl) 



R2 



P.Q - PQ. ^ P{PyQ - PQy) 



g3 
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Let a, Q{x,a) = 0, be a pole of Ri{x,y) of order ordQ(-Ri) = ordQ,(P/Q) = 
-ord„(g) = / > 1. We have ord„((P^Q-PQ^)g-2) < 21 and ord„(P(Pj,g- 
PQy)Q''^) = Sl + oidaiPyQ-PQy) = 2Z + 1 since ord^(P) = 0, ord„(Pj,Q) < 
—I, and ordQ,(PQy) = —/ + !. So orda(P2) =21 + 1. In general, the same 
argument shows that ordQ,(Pfc+i) = 2 ■ ordQ,(Pfc) + 1. Hence ordQ(Pfc) = 
2^~^l + 2^^^ — 1, k = 1,2, . . .. This is a strictly increasing sequence and we 
conclude again, since Rq, Ri, R2, ■ ■ ■ are hnearly independent over C(x), that 
A is not D-finite. ■ 

Proposition |21 and Theorem El show that I{x) is not D-finite and we can 
summarize the results of this section in the following corollary. 

Corollary 4 The sequences (in), (Com„), and (s„) are not P-recursive. 



3 Asymptotics 

We turn now to the computation of an asymptotic expansion for the numbers 
Sn- We will prove that: 

Theorem 5 

s„ = - 1 - - + + O . 

e^ \ n n[n — 1) I 



Our methods are such that, in principle, higher order terms could be obtained 
as a matter of brute force computation. In order to carry out this expansion 
we will first consider permutations which may not be simple, but whose non- 
trivial blocks all have length greater than some fixed value m. We will apply 
inclusion-exclusion arguments (dressed in the form of generating functions 
in]), an argument which allows us to reduce the number of terms considered, 
and a bootstrapping approach. 

The case m = 2, was already considered by Kaplansky 7]. Permutations 
of this type are those in which no two elements consecutive in position are 
also consecutive in value (in either order). These were called irreducible 
permutations by Atkinson and Stitt 0, but there is no standard terminology 
in the field. Indeed the permutations that we have referred to as plus and 
minus indecomposable have also been called irreducible in other contexts. 
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An amusing equivalent form for the case m = 2 is that the number of such 
permutations is also the number of ways of placing n mutually non-attacking 
krooks on an n X n chessboard. A krook is a piece which can move either like 
a king, or a rook in chess. Kaplansky's expansion is: 



In fact he derives asymptotic forms for the number of permutations contain- 
ing exactly r blocks of length 2 for any r. Our methods parallel his, and 
could also be used to derive such detailed information. 

The decomposition provided by Theorem^of a permutation into its maximal 
proper blocks represents a top down view of how non-simple permutations 
are constructed from simple ones. There is a corresponding bottom-up view 
that focuses on minimal blocks, put together in an arbitrary order. By a 
minimal block in vr we mean a non-singleton block in vr minimal with respect 
to inclusion. Note that the pattern of each minimal block is that of a simple 
permutation. Any permutation can be decomposed into minimal blocks and 
singletons, e.g., 3524716 = (3524) (7) (1) (6). However, this decomposition 
is not unique, for two essentially different reasons. The first one is that 
decompositions it = a[ai,a2, ■ ■ ■ , ak], where a is arbitrary and are simple, 
are not unique because it may be possible to coalesce singletons into simple 
blocks, or vice versa. Thus besides 3524716 = 2413[2413, 1, 1, 1] we also have 
3524716 = 3524716[1, 1, 1, 1, 1, 1, 1]. The second problem is that we require 
any two minimal blocks to be disjoint. While this is necessarily true whenever 
either of them has length more than 2, two minimal blocks of length 2 may 
intersect, as in 123. Thus we consider decompositions vr = a[ai, 02, • • • , 
where a is arbitrary and each a, is either 1, a simple permutation of length 
at least 4, or the identity permutation of length at least 2 or its reverse. We 
refer to blocks of the latter type as clusters in vr. 

By using clusters we have solved the second problem but the non-uniqueness 
remains and, moreover, we have introduced another source of it: consecutive 
(reversed) identical permutations may coalesce into longer (reversed) identi- 
cal permutations, as in 345612 = 21[1234, 12] = 231[12, 12, 12]. To remedy 
the non-uniqueness we introduce the notion of marking a permutation. A 
marked permutation (vr, M) consists of a permutation vr and a collection M 
of minimal blocks of it. A marked cluster in (tt, M) is a maximal chain of 
marked overlapping minimal blocks of length 2 (a marked cluster may be 
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a proper subset of a maximal cluster). Let Bi denote the set of all sim- 
ple permutations of length at least 4 and B2 denote the set of all identical 
permutations of length at least 2 and their reversals. Marking makes our 
decomposition unique: 

Theorem 6 Let X be the set of all marked permutations (tt, M) and Y be the 

set of all sequences {a;ai,a2,..., ak) where a is any permutation of length 
k > 1 and ai G {1} U ,Bi U i32. There is a bijection between the sets X and Y 
such that if (vr, M) (cr; ai, a2, ■ ■ ■ , CYk), where r of the ai belong to B\ and 
s of them to B2, then 

TT = a[ai, a2,..., ak] 
and \M\ = r + I — s where I is the total length of the a^ belonging to B2. 

Proof: Given a marked permutation, collapse its marked minimal blocks 
of length at least 4 and its marked clusters into singletons. This gives the 
permutation a. If the i-th term of a was not obtained by collapse then 
ttj = 1, otherwise a^ equals to the corresponding element of BiVJ B2- Since 
each ai G Bi contributes 1 to \M\ and each a^ G B2 of length m contributes 
m — 1, we have \M\ — r + I — s. It is clear that tt = crfai, a2, ■ ■ ■ , ak] and 
that (tt, M) can be uniquely recovered from {a;ai,a2,---, ak)- u 

Now suppose m to be some fixed value (we will later make choices of m 
suitable for our purposes, but will always assume that m > 2 since smaller 
values of m are trivial). Each permutation tt has an associated collection 
Bm{7r) consisting of the minimal blocks of tt whose length is less than or 
equal to m. So, if tt is simple and of length greater than m, Bm{T[) is empty, 
while for tt = 5672413, ^2(77) = {56, 67}, and ^4(77) = {56, 67, 2413}. An 
m-marking of tt is simply a subset of B^{7r). We consider the generating 
function: 

Then of course F^{x, —1) is the ordinary generating function for permuta- 
tions all of whose non-singleton blocks have length greater than m. 

We remark that -F„(x, i — 1) is the generating function where the coefficient 
of x'^t^ is precisely the number of permutations of length n with k minimal 
blocks of length less than or equal to m. 
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Let 



i=4 

We apply the bijection of Theorem El to marked permutations which contain 
no marked minimal blocks of length more than m. It follows that the gener- 
ating function of the corresponding permutations a G {1} Ui3i Ui32, in which 
X counts the length and v the contribution to |M|, is 



VX 

So: 

Fm{x, V) = y\k\ X + h vSm{x) 

from which it follows that: 



fc 



:= -1) = Y^k\[x- - Sm{x) \ . (4) 

Before using this equation to derive asymptotic information about s„ we 
digress briefly to show how it can be used to obtain an alternative derivation 
of (|21). Instead of using Sm[,x) in use S{x). This gives us foD^^^^i an 
ordinary generating function for permutations having no minimal block. The 
only such permutation is 1 so foo{x) = x. That is: 

2x^ 

X = F{x — S{x)) 

which yields (0) after applying to both sides. 

Now recall that fm{x) is the generating function for permutations all of whose 
blocks have length greater than m. In order to make use of these generating 
functions in the asymptotic computation of s„ we must determine a suitable 
value of m so that fm provides useful information about To that end the 
following lemma is useful. 

Lemma 7 // pn,k denotes the number of permutations of length n which 
contain a minimal block of length k then for any fixed positive integer c: 



— = 0{n 

k=c+2 "'■ 
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Proof: First observe that 

Pn,k < Sk{n - k + l){n - k + l)\ 

since the right hand side counts the number of ways to choose the structure 
of a block of length k, to choose its minimal element, and to arrange it with 
other elements, so it overcounts permutations with more than one such block. 

The estimate given then follows directly by using the fact that < /c!. Only 
the two extreme terms in the sum can have magnitude as large as 0{n^'^), 
and the remaining terms have magnitude 0{n~'^~^). Since there are fewer 
than n terms, the result follows. ■ 

So when seeking an asymptotic expansion of Sn/n\ with an error term of 
0(n~^~^) we may count instead the permutations which contain no blocks 
of length less than or equal to c + 2, or greater than or equal to n — c. 
In particular, as a direct consequence of the result quoted above due to 
Kaplansky [7j we obtain: 

Observation 8 



An alternative proof of this result follows from a more general theorem of 
Bender and Richmond |3j which provides the first order asymptotics of a 
class of series which include the inverse series of F{x). 

We will set as our goal to obtain the asymptotics of s„/n! with error term 
0{n^^). However, the technique we use is completely general, and could be 
applied, at the expense of a great deal of tedious computation, to any fixed 
error bound of this type. By the remarks above, we may ignore minimal 
block sizes between 5 and n — 3 inclusive. We first consider /4(x) which 
enumerates permutations having no minimal blocks of size less than or equal 
to 4. Recall that: 



•n 



-,+0{n-'). 




So: 



1 

n\ 
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-| oo 

= -yA;![r-'= 

1 " 
1 " 

= AE(-- 

Consider now any fixed value of / in equation (0). In order to obtain terms 
wliose order in n is or more, we need only consider the values 1 — 1 < i < I. 
Despite the fact that we sum over values of I running from through n, we 
may safely ignore the other terms. As we shall see in computing the three 
significant terms the summation over I does not affect the order of the terms. 

So, the three terms that we need to consider are: 





-I + 1)) + (6) 
(H + 2)H+l)/2 + /-2)). 

Each of these terms will be converted to the form: 

— — (an asymptotic expansion m n) . 

Since the first two and the first part of the third, are the same as those arising 
in the m = 2 case, we can make use of their known form, that is, use the 
asymptotics from Kaplansky's result, leaving only the term 

{n - /)!(-2)'-2(/ -2)fn- A (-2)' //(/ - !)(/ - 2) {n - l)\{n - 



n\ V^-2y l\ \ 4 n\{n-2l + 2)\ 



2y //(/-!)(/- 2) 



0{n 



l\ \ An{n - 1) 

Summing this expression over I gives —2e~'^/n{n — 1) + 0(r;,~^). 

Now we combine this additional term with Kaplansky's results to give the 
asymptotic expansion of [t"]/4(t) through three terms as: 

r]m = ^(i-—^,+o{n--')]. 

e"^ \ n{n — 1) / 
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Finally we use this in establishing the second order asymptotics of From 
Observation IHl applied to s„_i we obtain: 

sn-i = ^ (- + 0{n-')) . 

Furthermore, the number of permutations of length n containing a simple 
block of length n — 1 is precisely 4s„_i. Since, in computing the 1/n term in 
the expansion of s„ we can ignore contributions arising from blocks of length 
n — 2, and since the events of having a simple block of length from 2 to 4, 
and having a simple block of length (n — 1) are disjoint: 

Sn = [nMt)-4sn-i + 0{n-^n\y, 




We apply this bootstrap approach once more to get the second order be- 
haviour. We now know that: 

n n[n — 1) J 

\n[n — 1) j 

Furthermore there are 18s„_2 permutations of length n containing a simple 
block of length n — 2. However, of these 8s„_2 also contain a simple block of 
length 2. So: 

Sn = n/4(t)-4s„_i-10s„„2 + O(r2"3d) 

n\ f 4 2 ,A 

= - 1-- + ^ - + 0{n-')), 

e^ y n n[n — 1) J 

as we claimed at the beginning of this section. 

Finally, in this section we note that the asymptotic estimate of s„ is, as might 
be expected, a poor approximation. For example, S20 = 264111424634864638 
and our asymptotic estimate has a relative error of about 3.89 x 10~^. 

4 Congruences 

In this section we derive congruence properties of the numbers Com„ for 
the moduli 2*^ and 3 (from which follow similar congruences for Our 




15 



main tool is the following result that follows immediately from the Lagrange 
inversion formula. 



Lemma 9 



n ■ Com„ = [x"-^] I Y,{-lf{2\x + 3! 

vfc>0 



X 



For a prime p, let ordp(?2) denote the largest integer m such that divides 
n. As the following table shows, ord2(Com„) is unexpectedly large: 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


ord2(Com„) 





1 


1 


2 


2 


4 


4 


4 


4 


5 


5 


15 


13 


12 


12 



16 


17 


18 


19 


20 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


8 


8 


9 


9 


10 


10 


12 


12 


14 


14 


15 


15 


17 


17 


22 



In Theorem ^2 we give a lower bound on ord2(Com„) which is tight for 
infinitely many n and we completely characterize the values of n for which 
the equality is attained. 

For convenience we note the following result that follows directly from the 
well-known formula 



ordp(m! 



m 




m 






+ 




+ 


.P. 




p2 





Lemma 10 For all m, ord2{{m + l)\) > ^ where equality holds if and only 
if m = 1 or 2. Also, ord3(m!) < m — 1 for all m. 



Theorem 11 Let m = [n/2\. Then 

ord2(Com„) > 



n — 1 



Equality holds if and only if (^^^ is odd and this happens if and only if the 
binary expansion of m has no two consecutive unit digits. 
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Proof: Let the numbers bk, k > 0, he defined by 

k>0 k>0 

Thus bo = 1 and for A; > 1, 

bk= E (-l)^-(ci + l)!-(c2 + l)!-...-(c, + l)!. 

Cl,C2,...,Cs>l 

C1+C2H hCs=fc 

By Lemma ini 

n ■ Com„ = bk^bk^ ■■■bk„- 

fci,fc2,...,fc„>0 
ki+k2A \-k„=n—l 

By Lemma ITUl ord2((c + 1)!) > c/2 for all c. Hence, for all k and n, 

ord2(fefe) > - and ord2(r2 ■ Com„) > — - — . 

In particular, for odd n we have ord2(Com„) = ord2(n ■ Com^) > (n — l)/2. 

To obtain the more exact result of the theorem we need the following better 
estimates for ord2(6A:): 

{= k/2 for even k; 

= {k + l)/2 for A; = 1 mod 4; 
> (A; + l)/2 for A; = 3 mod 4. 

To prove them we look more closely at the sum for bk- Suppose first that 
k is even. Then the sum has exactly one summand with ord2 equal to k/2, 
namely that with ci = C2 = . . . = Ck/2 = 2 (by LemmalTUl ord2((c+l)!) = c/2 
only if c = 2), and the other summands have ord2 bigger than k/2. Hence 
ord2(&fc) = k/2. Now suppose that k is odd. Then each summand has an odd 
number of odd q's. The summands t with three and more odd q's satisfy 
ord2(t) > (A; + 3)/2 (each odd q contributes 1/2 to A;/2). The same is true if t 
has only one odd Cj but that q is not 1 (by LemmafTUl ord2((c+l)!) > (c+3)/2 
for odd c > 1), or if some even q is not 2 fLemma llUp. The remaining 
summands t, in which q = 2 with multiplicity {k — l)/2 and once q = 1, 
satisfy ord2(t) = {k + l)/2. We see that, for odd k, ord2(6fc) = {k + l)/2 
if and only if the number of the remaining summands is odd. This number 
equals (A;-l)/2 + l = (A; + l)/2. So ord2(6fc) = (A; + l)/2 if and only if A; = 1 
mod 4. 
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Let n = 2m + 1 be odd. If s is a summand of the above sum for n ■ Com„, 
then ord2(s) = (n — l)/2 if and only if aU ki in s are even; other summands 
t have ord2(t) > (n — l)/2. It follows that ord2(Com„) = (n — l)/2 if and 
only if the number of the former summands s is odd. This number equals 




, m J 



Let n = 2m be even. We know that ord2(6A;) = k/2 for even k and ord2(6fe) > 
(A; + l)/2 for odd k. In the sum for n ■ Com„, every composition fci + A;2 + 
■■■-\-kn — n — loin — 1 has an odd number of odd parts. For any i-tuple 
Zi, I2, ■ ■ ■ , h, where t and all li are odd and li + ■ ■ ■ + It < n — 1, we let 
S{li, I2, ■ ■ ■ , h) denote the sum of those bk^bk^ with ki + k2 + ■ ■ ■ + kn — 

n — 1 in which k^ — l^, 1 < i < t, and k^ is even for i > t. It follows that 

n ■ Com„ = J2 Q'S'(/i,/2, ■■■,k) 

where we sum over all mentioned i-tuples /i, /2, • • • , k- By the properties of 
ord2 and of the numbers bk, ord2(S'(/i, h, ■ ■ ■ , k)) > {n + t — l)/2. Also, for 
odd t we have ord2((^)) = ordsCf = ord2(n) - ord2(i) +ord2(('^rJ)) > 

ord2(n), and ord2((")) = ord2(n). It follows that ord2(Com„) > n/2 and, 
moreover, ord2(Com„) — n/2 if and only if 

ord2 I =n/2. 

\l<n, I odd / 

In the last sum still many summands have ord2 bigger than n/2: if / = 3 
mod 4 then ord(iS'(/)) > n/2. On the other hand, if / = 1 mod 4 then each 
summand bibk2 ■ ■ ■ b^^ in S{1) has OTd2{bibk2 . . . 6^^) = n/2. We conclude that 
ord2(Com„) = n/2 if and only if the number c(n) of compositions of n — 1 
into n parts, where the first part is = 1 mod 4 and the remaining n — 1 parts 
are even (zero parts are allowed), is odd. We have 



cm 



[a^""']:r^ • TT^-^ = K"']7; ^^wT mod 2 
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3m — l\ _ 3m / 3m — 1 
— 1 j m \m — 1 ^ 
'3m\ 
.ml' 



mod 2 



It was noted by Kummer see also Singmaster jTTj, that oidp{{^'^^^) is 
equal to the number of carries required when adding a and h in the p-ary 
notation. Applying this for p = 2, a = m, and b = 2m, we get the stated 
criterion. ■ 

Corollary 12 For all n > 3, 

^ f 2 mod 2("-i)/2 for odd n; 
~ 1 —2 mod 2"/^ /or etJen n. 

Let 

n + 1 \ n y 

be the nth Catalan number. 

Proposition 13 For all n, Com„ = C„_i mod 3. 

Proof: We have, for every non-negative integer /c, 

(2!x + 3!x2 + ■ ■ -f = (2a;)^ + 3ak{x) 
with ak{x) G Z[[,t]]. Thus 

E(-l)'(2!a: + 3!x^ + ■■■)'= = + 3 

fc>0 i + 

^ 3b{x) 



l + 2x 



with e Z[[a;]]. Let m = ord3(?T,). Since ord3(/c!) < /c — 1 for every A; 
fLemmalTUj). we have 

ord3(3^Q) > m + 1 for /c = l,2,...,n. 
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By Lemma ini 



n ■ Com„ = [x""'] — + 3b{x 

\ ±. I iLiX 



1 \" 1 

„n-iW , Q^/^^\ = Fa-n-ii^ -mod 3™+^ 



(l + 2x) 

-2 



n-i/2n-2^ 



Y n — 1 y 

Canceling in the last congruence the common factor 3™, we get 

^ . Com.. . <^P« - - 2) mod 3. 

3™ 3"^ V ^ - 1 / 3™ V ^ - 1 y 

Since n/3™ ^ mod 3, we can divide by it and get 

1 f2n- 2\ 
Com„ = — mod 3. 

n \ n — 1 



Corollary 14 For all n > 2, 

Sn = -Cn-i + (-1)" mod 3. 



5 Concluding remarks 

The simplicity property for permutations does not seem to have been stud- 
ied until very recently jTUl HI- We have begun the study of the numbers 
Sn by showing that they are not P-recursive, giving the first few terms of 
their asymptotic expansion, and showing that they satisfy some unexpected 
congruence properties. 

These results suggest a number of natural continuations. Although, in prin- 
ciple, we could obtain more terms of the asymptotic expansion the entire 
expansion remains elusive, and computing it seems to be rather a difficult 
problem. On the other hand we have some computational evidence to suggest 
that the sequence Com^ has additional congruence properties, particularly 
with respect to odd primes. 

We suggest also some algorithmic problems that are natural counterparts to 
the enumerative results: 
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• How can one efficiently generate simple permutations in lexicographic 
order? 

• Is it possible to generate simple permutations uniformly at random in 
worst-case linear time per permutation? 

• How efficiently can one recognise a simple permutation? 

With regards to the final question, there is a natural dynamic programming 
algorithm that achieves the task in Oln"^) time; so the issue is whether one 
can do better. 
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